Subgenome-aware analyses reveal the genomic consequences of ancient allopolyploid hybridizations throughout the cotton family.

Malvaceae comprise some 4,225 species in 243 genera and nine subfamilies and include economically important species, such as cacao, cotton, durian, and jute, with cotton an important model system for studying the domestication of polyploids. Here, we use chromosome-level genome assemblies from representatives of five or six subfamilies (depending on the placement of Ochroma) to differentiate coexisting subgenomes and their evolution during the family's deep history. The results reveal that the allohexaploid Helicteroideae partially derive from an allotetraploid Sterculioideae and also form a component of the allodecaploid Bombacoideae and Malvoideae. The ancestral Malvaceae karyotype consists of 11 protochromosomes. Four subfamilies share a unique reciprocal chromosome translocation, and two other subfamilies share a chromosome fusion. DNA alignments of single-copy nuclear genes do not yield the same relationships as inferred from chromosome structural traits, probably because of genes originating from different ancestral subgenomes. These results illustrate how chromosome-structural data can unravel the evolutionary history of groups with ancient hybrid genomes.

A sequence-aware merger of genomic structural variations at population scale.

Merging structural variations (SVs) at the population level presents a significant challenge, yet it is essential for conducting comprehensive genotypic analyses, especially in the era of pangenomics. Here, we introduce PanPop, a tool that utilizes an advanced sequence-aware SV merging algorithm to efficiently merge SVs of various types. We demonstrate that PanPop can merge and optimize the majority of multiallelic SVs into informative biallelic variants. We show its superior precision and lower rates of missing data compared to alternative software solutions. Our approach not only enables the filtering of SVs by leveraging multiple SV callers for enhanced accuracy but also facilitates the accurate merging of large-scale population SVs. These capabilities of PanPop will help to accelerate future SV-related studies.

Microbial and metabolic profiles unveil mutualistic microbe-microbe interaction in obesity-related colorectal cancer.

Obesity is a risk factor for colorectal cancer (CRC), and the involvement of gut microbiota in the pathogenesis of obesity and CRC is widely recognized. However, the landscape of fecal microbiome and metabolome distinguishing patients with obesity-related CRC from obesity remains unknown. Here, we utilize metagenomic sequencing and metabolomics from 522 patients with CRC and healthy controls to identify the characteristics of obese CRC. Our integrated analysis reveals that obesity-related CRC is characterized by elevated Peptostreptococcus stomatis, dysregulated fatty acids and phospholipids, and altered Kyoto Encyclopedia of Genes and Genomes pathways involving glycerophospholipid metabolism and lipopolysaccharide synthesis. Correlation analysis unveils microbial interactions in obesity, where the probiotic Faecalibacterium prausnitzii and the tumor-promoting species P. stomatis may engage in cross-feeding, thereby promoting tumorigenesis. In vitro experiments affirm enhanced growth under cross-feeding conditions. The mutualistic microbe-microbe interaction may contribute to the association between obesity and elevated CRC risk. Additionally, diagnostic models incorporating BMI-specific microbial biomarkers display promise for precise CRC screening.

Chromosome-scale genome assembly of Lepus oiostolus (Lepus, Leporidae).

Lepus oiostolus (L. oiostolus) is a species endemic to the Qinghai-Tibet Plateau. However, the absence of a reference genome limits genetic studies. Here, we reported a high-quality L. oiostolus genome assembly, with scaffolds anchored to 24 chromosomes and a total assembled length of 2.80 Gb (contig N50 = 64.25 Mb). Genomic annotation uncovered 22,295 protein-coding genes and identified 49.84% of the sequences as transposable elements. Long interspersed nuclear elements (LINEs) constitute a high proportion of the genome. Our study is at the first time to report the chromosome-scale genome for the species of the L. oiostolus. It provides a valuable genomic resource for future research on the evolution of the Leporidae.

Autophagy protein ATG-16.2 and its WD40 domain mediate the beneficial effects of inhibiting early-acting autophagy genes in C. elegans neurons.

While autophagy genes are required for lifespan of long-lived animals, their tissue-specific roles in aging remain unclear. Here, we inhibited autophagy genes in Caenorhabditis elegans neurons, and found that knockdown of early-acting autophagy genes, except atg-16.2, increased lifespan, and decreased neuronal PolyQ aggregates, independently of autophagosomal degradation. Neurons can secrete protein aggregates via vesicles called exophers. Inhibiting neuronal early-acting autophagy genes, except atg-16.2, increased exopher formation and exopher events extended lifespan, suggesting exophers promote organismal fitness. Lifespan extension, reduction in PolyQ aggregates and increase in exophers were absent in atg-16.2 null mutants, and restored by full-length ATG-16.2 expression in neurons, but not by ATG-16.2 lacking its WD40 domain, which mediates noncanonical functions in mammalian systems. We discovered a neuronal role for C. elegans ATG-16.2 and its WD40 domain in lifespan, proteostasis and exopher biogenesis. Our findings suggest noncanonical functions for select autophagy genes in both exopher formation and in aging.

Multiple independent losses of the biosynthetic pathway for two tropane alkaloids in the Solanaceae family.

Hyoscyamine and scopolamine (HS), two valuable tropane alkaloids of significant medicinal importance, are found in multiple distantly related lineages within the Solanaceae family. Here we sequence the genomes of three representative species that produce HS from these lineages, and one species that does not produce HS. Our analysis reveals a shared biosynthetic pathway responsible for HS production in the three HS-producing species. We observe a high level of gene collinearity related to HS synthesis across the family in both types of species. By introducing gain-of-function and loss-of-function mutations at key sites, we confirm the reduced/lost or re-activated functions of critical genes involved in HS synthesis in both types of species, respectively. These findings indicate independent and repeated losses of the HS biosynthesis pathway since its origin in the ancestral lineage. Our results hold promise for potential future applications in the artificial engineering of HS biosynthesis in Solanaceae crops.

KIF15 knockdown inhibits colorectal cancer proliferation and migration through affecting the ubiquitination modification of NRAS.

As one of the most common malignancies, colorectal cancer (CRC) requires a thorough understanding of the mechanisms that promote its development and the discovery of new therapeutic targets. In this study, immunohistochemical staining confirmed significantly higher expression levels of KIF15 in CRC. qPCR and western blot results demonstrated the effective suppression of KIF15 mRNA and protein expression by shKIF15. Downregulation of KIF15 inhibited the proliferation and migration of CRC cells while promoting apoptosis. In addition, evidence from the xenograft experiments in nude mice demonstrated that KIF15 knockdown also suppressed tumor growth. Through bioinformatics analysis, the downstream molecular NRAS and Rac signaling pathway associated with KIF15 were identified. KIF15 knockdown was found to inhibit NRAS expression and disrupt Rac signaling pathway. Moreover, WB and Co-IP assays revealed that KIF15 reduced the ubiquitination modification of NRAS protein by interacting with the E3 ligase MDM2, thereby enhancing NRAS protein stability. Functionally, NRAS knockdown was shown to inhibit cell proliferation and migration. In conclusion, KIF15 promoted CRC progression by regulating NRAS expression and Rac signaling pathway.

Subgenomic Stability of Progenitor Genomes During Repeated Allotetraploid Origins of the Same Grass Brachypodium hybridum.

Both homeologous exchanges and homeologous expression bias are generally found in most allopolyploid species. Whether homeologous exchanges and homeologous expression bias differ between repeated allopolyploid speciation events from the same progenitor species remains unknown. Here, we detected a third independent and recent allotetraploid origin for the model grass Brachypodium hybridum. Our homeologous exchange with replacement analyses indicated the absence of significant homeologous exchanges in any of the three types of wild allotetraploids, supporting the integrity of their progenitor subgenomes and the immediate creation of the amphidiploids. Further homeologous expression bias tests did not uncover significant subgenomic dominance in different tissues and conditions of the allotetraploids. This suggests a balanced expression of homeologs under similar or dissimilar ecological conditions in their natural habitats. We observed that the density of transposons around genes was not associated with the initial establishment of subgenome dominance; rather, this feature is inherited from the progenitor genome. We found that drought response genes were highly induced in the two subgenomes, likely contributing to the local adaptation of this species to arid habitats in the third allotetraploid event. These findings provide evidence for the consistency of subgenomic stability of parental genomes across multiple allopolyploidization events that led to the same species at different periods. Our study emphasizes the importance of selecting closely related progenitor species genomes to accurately assess homeologous exchange with replacement in allopolyploids, thereby avoiding the detection of false homeologous exchanges when using less related progenitor species genomes.

Scattered differentiation of unlinked loci across the genome underlines ecological divergence of the selfing grass Brachypodium stacei.

Ecological divergence without geographic isolation, as an early speciation process that may lead finally to reproductive isolation through natural selection, remains a captivating topic in evolutionary biology. However, the pattern of genetic divergence underlying this process across the genome may vary between species and mating systems. Here, we present evidence that Brachypodium stacei, an annual and highly selfing grass model species, has undergone sympatric ecological divergence without geographic isolation. Genomic, transcriptomic, and metabolomic analyses together with lab experiments mimicking the two opposite environmental conditions suggest that diploid B. stacei populations have diverged sympatrically in two slopes characterized by distinct biomes at Evolution Canyon I (ECI), Mount Carmel, Israel. Despite ongoing gene flow, primarily facilitated by seed dispersal, the level of gene flow has progressively decreased over time. This local adaptation involves the scattered divergence of many unlinked loci across the total genome that include both coding genes and noncoding regions. Additionally, we have identified significant differential expressions of genes related to the ABA signaling pathway and contrasting metabolome composition between the arid- vs. forest-adapted B. stacei populations in ECI. These results suggest that multiple small loci involved in environmental responses act additively to account for ecological adaptations by this selfing species in contrasting environments.

Proteostasis is differentially modulated by inhibition of translation initiation or elongation.

Recent work has revealed an increasingly important role for mRNA translation in maintaining proteostasis. Here, we use chemical inhibitors targeting discrete steps of translation to compare how lowering the concentration of all or only translation initiation-dependent proteins rescues Caenorhabditis elegans from proteotoxic stress. We systematically challenge proteostasis and show that pharmacologically inhibiting translation initiation or elongation elicits a distinct protective profile. Inhibiting elongation protects from heat and proteasome dysfunction independently from HSF-1 but does not protect from age-associated protein aggregation. Conversely, inhibition of initiation protects from heat and age-associated protein aggregation and increases lifespan, dependent on hsf-1, but does not protect from proteotoxicity caused by proteasome dysfunction. Surprisingly, we find that the ability of the translation initiation machinery to control the concentration of newly synthesized proteins depends on HSF-1. Inhibition of translation initiation in wild-type animals reduces the concentration of newly synthesized proteins but increases it in hsf-1 mutants. Our findings suggest that the HSF-1 pathway is not only a downstream target of translation but also directly cooperates with the translation initiation machinery to control the concentration of newly synthesized proteins to restore proteostasis.

Function and mechanism of MCM8 in the development and progression of colorectal cancer.

Colorectal cancer (CRC) has become a global health problem which has almost highest morbidity and mortality in all types of cancers. This study aimed to uncover the biological functions and underlying mechanism of MCM8 in the development and progression of CRC. The expression level of MCM8 was found to be upregulated in CRC tissues and significantly associated with tumor grade and patients' survival. Knocking down MCM8 expression in CRC cells could restrain cell growth and cell motility while promoting cell apoptosis in vitro, as well as inhibit tumor growth in xenograft mice model. Based on the RNA screening performing on CRC cells with or without MCM8 knockdown and the following IPA analysis, CHSY1 was identified as a potential target of MCM8 in CRC, whose expression was also found to be higher in tumor tissues than in normal tissues. Moreover, it was demonstrated that MCM8 may regulate the expression of CHSY1 through affecting its NEDD4-mediated ubiquitination, both of which synergistically execute tumor promotion effects on CRC. In conclusion, the outcomes of our study showed the first evidence that MCM8 act as a tumor promotor in CRC, and may be a promising therapeutic target of CRC treatment.

Multi-omics data provide insight into the adaptation of the glasshouse plant Rheum nobile to the alpine subnival zone.

Subnival glasshouse plants provide a text-book example of high-altitude adaptation with reproductive organs enclosed in specialized semi-translucent bracts, monocarpic reproduction and continuous survival under stress. Here, we present genomic, transcriptomic and metabolomic analyses for one such plant, the Noble rhubarb (Rheum nobile). Comparative genomic analyses show that an expanded number of genes and retained genes from two recent whole-genome duplication events are both relevant to subnival adaptation of this species. Most photosynthesis genes are downregulated within bracts compared to within leaves, and indeed bracts exhibit a sharp reduction in photosynthetic pigments, indicating that the bracts no longer perform photosynthesis. Contrastingly, genes related to flavonol synthesis are upregulated, providing enhanced defense against UV irradiation damage. Additionally, anatomically abnormal mesophyll combined with the downregulation of genes related to mesophyll differentiation in bracts illustrates the innovation and specification of the glass-like bracts. We further detect substantial accumulation of antifreeze proteins (e.g. AFPs, LEAs) and various metabolites (e.g. Proline, Protective sugars, procyanidins) in over-wintering roots. These findings provide new insights into subnival adaptation and the evolution of glasshouse alpine plants.

Evolutionary origin of genomic structural variations in domestic yaks.

Yak has been subject to natural selection, human domestication and interspecific introgression during its evolution. However, genetic variants favored by each of these processes have not been distinguished previously. We constructed a graph-genome for 47 genomes of 7 cross-fertile bovine species. This allowed detection of 57,432 high-resolution structural variants (SVs) within and across the species, which were genotyped in 386 individuals. We distinguished the evolutionary origins of diverse SVs in domestic yaks by phylogenetic analyses. We further identified 334 genes overlapping with SVs in domestic yaks that bore potential signals of selection from wild yaks, plus an additional 686 genes introgressed from cattle. Nearly 90% of the domestic yaks were introgressed by cattle. Introgression of an SV spanning the KIT gene triggered the breeding of white domestic yaks. We validated a significant association of the selected stratified SVs with gene expression, which contributes to phenotypic variations. Our results highlight that SVs of different origins contribute to the phenotypic diversity of domestic yaks.

A chromosome-scale Rhubarb (Rheum tanguticum) genome assembly provides insights into the evolution of anthraquinone biosynthesis.

Rhubarb is the collective name for various perennial plants from the genus Rheum L. and the Polygonaceae family. They are one of the most ancient, commonly used, and important herbs in traditional Chinese medicine. Rhubarb is a major source of anthraquinones, but how they are synthesized remains largely unknown. Here, we generate a genome sequence assembly of one important medicinal rhubarb R. tanguticum at the chromosome level, with 2.76 Gb assembled into 11 chromosomes. The genome is shaped by two recent whole-genome duplication events and recent bursts of retrotransposons. Metabolic analyses show that the major anthraquinones are mainly synthesized in its roots. Transcriptomic analysis reveals a co-expression module with a high correlation to anthraquinone biosynthesis that includes key chalcone synthase genes. One CHS, four CYP450 and two BGL genes involved in secondary metabolism show significantly upregulated expression levels in roots compared with other tissues and clustered in the co-expression module, which implies that they may also act as candidate genes for anthraquinone biosynthesis. This study provides valuable insights into the genetic bases of anthraquinone biosynthesis that will facilitate improved breeding practices and agronomic properties for rhubarb in the future.

Depicting the landscape of gut microbial-metabolic interaction and microbial-host immune heterogeneity in deficient and proficient DNA mismatch repair colorectal cancers.

BACKGROUND: Accumulating evidence has indicated the role of gut microbiota in remodeling host immune signatures, but various interplays underlying colorectal cancers (CRC) with deficient DNA mismatch repair (dMMR) and proficient DNA mismatch repair (pMMR) remain poorly understood. This study aims to decipher the gut microbiome-host immune interactions between dMMR and pMMR CRC. METHOD: We performed metagenomic sequencing and metabolomic analysis of fecal samples from a cohort encompassing 455 participants, including 21 dMMR CRC, 207 pMMR CRC, and 227 healthy controls. Among them, 50 tumor samples collected from 5 dMMR CRC and 45 pMMR CRC were conducted bulk RNA sequencing. RESULTS: Pronounced microbiota and metabolic heterogeneity were identified with 211 dMMR-enriched species, such as Fusobacterium nucleatum and Akkermansia muciniphila, 2 dMMR-depleted species, such as Flavonifractor plautii, 13 dMMR-enriched metabolites, such as retinoic acid, and 77 dMMR-depleted metabolites, such as lactic acid, succinic acid, and 2,3-dihydroxyvaleric acid. F. plautii was enriched in pMMR CRC and it was positively associated with fatty acid degradation, which might account for the accumulation of dMMR-depleted metabolites classified as short chain organic acid (lactic acid, succinic acid, and 2,3-dihydroxyvaleric acid) in pMMR CRC. The microbial-metabolic association analysis revealed the characterization of pMMR CRC as the accumulation of lactate induced by the depletion of specific gut microbiota which was negatively associated with antitumor immune, whereas the nucleotide metabolism and peptide degradation mediated by dMMR-enriched species characterized dMMR CRC. MMR-specific metabolic landscapes were related to distinctive immune features, such as CD8+ T cells, dendritic cells and M2-like macrophages. CONCLUSIONS: Our mutiomics results delineate a heterogeneous landscape of microbiome-host immune interactions within dMMR and pMMR CRC from aspects of bacterial communities, metabolic features, and correlation with immunocyte compartment, which infers the underlying mechanism of heterogeneous immune responses.

The genome sequence and demographic history of Przewalskia tangutica (Solanaceae), an endangered alpine plant on the Qinghai-Tibet Plateau.

To adapt to high-altitude habitats, many alpine plants develop self-compatible breeding systems from outcrossing. The genetic bases for this shift and the resulting demographic consequences remain largely unexplored. Here, we present a high-quality, chromosome-level genome assembly of the monotypic and endangered alpine perennial Przewalskia tangutica (Solanaceae) occurring on the Qinghai-Tibet Plateau (QTP). Our assembled genome is approximately 3 Gb, with a contig N50 size of 17 Mb, and we identified one lineage-specific whole-genome duplication. We found that the gametophytic self-incompatibility (GSI) syntenic locus to the other obligate outcrossing Solanaceae species was broken by the inserted the long terminal repeats, and changes in the flower-specific expression of the homologous genes, and the linked GSI genes in this species. Such changes may have led to its self-compatibility. We identified three deeply diverged lineages in the central distribution of this species, and the gene flow between them was weak but continuous. All three lineages diverged and decreased their population sizes since the largest glaciations occurred in the QTP approximately 720-500 thousand years ago. In addition, we identified one obvious hybrid population between two lineages, suggesting that genetic exchanges between and within lineages still occur. Our results provide insights into evolutionary adaptation through facultative self-pollination and demographic consequences of this alpine rare species in arid habitats.

Genome-scale angiosperm phylogenies based on nuclear, plastome, and mitochondrial datasets.

Angiosperms dominate the Earth's ecosystems and provide most of the basic necessities for human life. The major angiosperm clades comprise 64 orders, as recognized by the APG IV classification. However, the phylogenetic relationships of angiosperms remain unclear, as phylogenetic trees with different topologies have been reconstructed depending on the sequence datasets utilized, from targeted genes to transcriptomes. Here, we used currently available de novo genome data to reconstruct the phylogenies of 366 angiosperm species from 241 genera belonging to 97 families across 43 of the 64 orders based on orthologous genes from the nuclear, plastid, and mitochondrial genomes of the same species with compatible datasets. The phylogenetic relationships were largely consistent with previously constructed phylogenies based on sequence variations in each genome type. However, there were major inconsistencies in the phylogenetic relationships of the five Mesangiospermae lineages when different genomes were examined. We discuss ways to address these inconsistencies, which could ultimately lead to the reconstruction of a comprehensive angiosperm tree of life. The angiosperm phylogenies presented here provide a basic framework for further updates and comparisons. These phylogenies can also be used as guides to examine the evolutionary trajectories among the three genome types during lineage radiation.

Turfgrass-dependent mycotrophic change enhances soil deterioration in dry, cold and high-alkali environments.

AIMS: Soil quality is undergoing severe degradation under anthropogenic effects. Different methods of land management have been implemented for soil reclamation, such as turfing. Although widely accepted to improve soil quality, turfing in specific environments may also culminate in soil deterioration. We aim to know how turfing impacts soils by changing mycobiomes. METHODS AND RESULTS: The soil physicochemical properties and ITS metabarcoding were used to investigate mycobiome diversity and eco-function differences between the eudicot Dianthus plumarius and the monocot Poa pratensis in dry, cold, and high-alkali soil. The effects of plantation and the rhizosphere (e.g. root exudates) were tested. We showed that the change in soil mycobiomes in different planted bulk soils and rhizospheres could mainly be attributed to species turnover, with minor nestedness. Unexpectedly, the soil deteriorates more following turfing. The increasing saprotrophs in planted bulk soil were more marked in the monocot than in the eudicot, even the rhizosphere effect alleviated saprotrophic risks in the rhizosphere. CONCLUSIONS: Turfing deteriorates the health of high-alkali soil by reducing nitrification, and upshift the soil saprotrophs in a dry and cold environment.

Integrated metagenomic and metabolomic analysis reveals distinct gut-microbiome-derived phenotypes in early-onset colorectal cancer.

OBJECTIVE: The incidence of early-onset colorectal cancer (EO-CRC) is steadily increasing. Here, we aimed to characterise the interactions between gut microbiome, metabolites and microbial enzymes in EO-CRC patients and evaluate their potential as non-invasive biomarkers for EO-CRC. DESIGN: We performed metagenomic and metabolomic analyses, identified multiomics markers and constructed CRC classifiers for the discovery cohort with 130 late-onset CRC (LO-CRC), 114 EO-CRC subjects and age-matched healthy controls (97 LO-Control and 100 EO-Control). An independent cohort of 38 LO-CRC, 24 EO-CRC, 22 LO-Controls and 24 EO-Controls was analysed to validate the results. RESULTS: Compared with controls, reduced alpha-diversity was apparent in both, LO-CRC and EO-CRC subjects. Although common variations existed, integrative analyses identified distinct microbiome-metabolome associations in LO-CRC and EO-CRC. Fusobacterium nucleatum enrichment and short-chain fatty acid depletion, including reduced microbial GABA biosynthesis and a shift in acetate/acetaldehyde metabolism towards acetyl-CoA production characterises LO-CRC. In comparison, multiomics signatures of EO-CRC tended to be associated with enriched Flavonifractor plauti and increased tryptophan, bile acid and choline metabolism. Notably, elevated red meat intake-related species, choline metabolites and KEGG orthology (KO) pldB and cbh gene axis may be potential tumour stimulators in EO-CRC. The predictive model based on metagenomic, metabolomic and KO gene markers achieved a powerful classification performance for distinguishing EO-CRC from controls. CONCLUSION: Our large-sample multiomics data suggest that altered microbiome-metabolome interplay helps explain the pathogenesis of EO-CRC and LO-CRC. The potential of microbiome-derived biomarkers as promising non-invasive tools could be used for the accurate detection and distinction of individuals with EO-CRC.

Multi-omics reveal differentiation and maintenance of dimorphic flowers in an alpine plant on the Qinghai-Tibet Plateau.

Dimorphic flowers growing on a single individual plant play a critical role in extreme adaption and reproductive assurance in plants and have high ecological and evolutionary significance. However, the omics bases underlying such a differentiation and maintenance remain largely unknown. We aimed to investigate this through genomic, transcriptome and metabolomic analyses of dimorphic flowers in an alpine biennial, Sinoswertia tetraptera (Gentianaceae). A high-quality chromosome-level genome sequence (903 Mb) was first assembled for S. tetraptera with 31,359 protein-coding genes annotated. Two rounds of recent independent whole-genome duplication (WGD) were revealed. Numerous genes from the recent species-specific WGD were found to be differentially expressed in the two types of flowers, and this may have helped contribute to the origin of this innovative trait. The genes with contrasting expressions between flowers were related to biosynthesis of hormones, floral pigments (carotenoids and flavonoids) and iridoid compounds, which are involved in both flower development and colour. Metabolomic analyses similarly suggested differential concentrations of these chemicals in the two types of flowers. The expression interactions between multiple genes may together lead to contrasting morphology and chemical concentration and open versus closed pollination of the dimorphic flowers in this species for reproductive assurance.